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Abstract 

O 

We consider a half-duplex wireless relay network with hybrid-automatic retransmission request 
(HARQ) and Rayleigh fading channels. In this paper, we analyze the outage probability of the multi- 
_ relay delay-limited HARQ system with opportunistic relaying scheme in decode-and-forward mode, in 

. which the best relay is selected to transmit the source's regenerated signal. A simple and distributed 

(N 

(N 
O 
(N 



relay selection strategy is proposed for multi-relay HARQ channels. Then, we utilize the non-orthogonal 

(N ■ 

f"*) ■ cooperative transmission between the source and selected relay for retransmitting of the source data 

toward the destination if needed, using space-time codes or beamforming techniques. We analyze the 
performance of the system. We first derive the cumulative density function (CDF) and probability 

■ density function (PDF) of the selected relay HARQ channels. Then, the CDF and PDF are used to 

■ determine the outage probability in the /-th round of HARQ. The outage probability is required to 
compute the throughput-delay performance of this half-duplex opportunistic relaying protocol. The 
packet delay constraint is represented by L, the maximum number of HARQ rounds. An outage is 
declared if the packet is unsuccessful after L HARQ rounds. Furthermore, closed-form upper-bounds 
on outage probability are derived and subsequently are used to investigate the diversity order of the 
system. Based on the derived upper-bound expressions, it is shown that the proposed schemes achieve 
the full spatial diversity order of + 1, where N is the number of potential relays. Our analytical 
results are confirmed by simulation results. 
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I. Introduction 

Cooperation among devices has been considered to provide diversity in wireless networks 
where fading may significantly affect single links [1]. Initial works have emphasized on relaying, 
where a cooperator node amplifies (or decodes) and forwards, possibly in a quantized fashion [2], 
the information from the source node in order to help decoding at the destination node [3], [4], 
[5]. The achieved throughput can be increased with the integration of cooperation and coding, 
i.e., by letting the cooperator send incremental redundancy to the destination [6]. In particular, 
it has been shown in [6] that coded cooperation achieves a diversity order of two, while decode- 
and-forward reaches only a diversity order one, when the transmissions of source and cooperator 
are orthogonal. The capacity of cooperative networks using both the decode-and-forward and 
coded cooperation has been extensively studied [6], [7] for simple networks with simple medium 
access control (MAC) protocols. In [8], a system with two transmission phases that makes use of 
convolutional codes is analyzed and characterized by means of partner choice and performance 
regions. Resource allocation for space-time coded cooperative networks has been studied in [9], 
[10], where the analysis of bit error rate and outage probability are also derived. Unfortunately, 
in cooperative relaying the diversity gain is increased at the expense of throughput loss due to the 
half-duplex constraint at relay nodes. Different methods have been proposed to recover this loss. 
In [11], successive relaying using repetition coding has been introduced for a two relay wireless 
network with flat fading. In [12], relay selection methods have been proposed for cooperative 
communication with decode-and-forward (DF) relaying. 

A prominent alternative to reducing the throughput loss in relay-aided transmission mecha- 
nisms is the combination of both ARQ and relaying. This approach would significantly reduce 
the half-duplex multiplexing loss by activating ARQ for rare erroneously decoded data packets, 
when they occur. Approaches targeting the joint design of ARQ and relaying in one common 
protocol have recently received more interest (see for instance [13], [14]). Motivated by the above 
suggestion, we investigate and analyze throughput efficient cooperative transmission techniques 
where both ARQ and relaying are jointly designed. To this end, a diversity effect can be 
introduced to a relay networks by simply allowing the nodes to maintain previously received 
information concerning each active message. Each time a message is retransmitted, either from 
a new node or from the same node, every node in the relay network will increase the amount 
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of resolution information it has about the message. Once a node has accumulated sufficient 
information it will be able to decode the message and can act as a relay and forward the 
message (as in decode-and-forward [3], [15]). This diversity effect can be viewed as a space- 
time generalization of the time-diversity effect of hybrid- automatic repeat request (HARQ) as 
described in [16]. Thus, the HARQ scheme which is used in this paper is a practical approach 
to designing wireless ad hoc networks that exploit the spatial diversity, which is achievable 
with relaying. The retransmitted packets could originate from any node that has overheard 
and successfully decoded the message. Current and future wireless networks based on packet 
switching use HARQ protocols at the link layer. Hence, the performance of HARQ protocols in 
relay channels has attracted recent research interest [17], [18], [13]. 

In this paper, we propose an efficient HARQ multi-relay protocol which leads to full spatial 
diversity. We assume a delay-limited network with the maximum number of HARQ rounds 
L, which represents the delay constraint. The protocol uses a form of incremental redundancy 
HARQ transmission with assistance from the selected relay via non-orthogonal transmission in 
the second transmission phase if the relay decodes the message before the destination. Note 
that by non-orthogonal transmission, we mean that the source node and the selected relay 
simultaneously retransmit the source data using space-time codes or beamforming techniques. 
We introduce a distributed relay selection scheme for HARQ multi-relay networks by using 
acknowledgment (ACK) or non- acknowledgment (NACK) signals transmitted by destination. 
Closed-form expressions are derived for the outage probability, defined as the probability of 
packet failure after L HARQ rounds, in half-duplex. For sufficiently high SNR, we derive a 
simple closed-form average outage probability expression for a HARQ system with multiple 
cooperating branches, and it is shown that the full diversity is achievable in the proposed HARQ 
relay networks. The simulations shows that the throughput of the relay channel is significantly 
larger than that of direct transmission for a wide range of signal-to-noise ratios (SNRs), target 
outage probabilities, delay constraints and relay numbers. 

The remainder of this paper is organized as follows: In Section II, the system model and 
protocol description are given. The performance analysis The closed-form expressions for the 
outage probability and asymptotic analysis of the system are presented in Section III, which are 
utilized for optimizing the system. In Section IV, the overall system performance is presented 
for different numbers of relays and channel conditions, and the correctness of the analytical 
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formulas are confirmed by simulation results. Conclusions are presented in Section V. 

Notations: The superscripts (■)*, (■)^, and (■)* stand for transposition, conjugate transposition, 
and element- wise conjugation, respectively. The expectation operation is denoted by IE{-}. The 
union and intersection of a collection of sets are denoted by |J and f], respectively. The symbol 
|x| is the absolute value of the scalar x, while [x]^ denotes max{a;, 0}. The logarithms log2 and 
log are the based two logarithm and the natural logarithm, respectively. 

IL System Model and Protocol Description 

Consider a network consisting of a source, one or more relays denoted i = 1,2, . . . , N, and 
one destination. The wireless relay network model is illustrated in Fig. [H It is assumed that 
each node is equipped with a single antenna. We consider symmetric channels and denote the 
source-to-destination, source-to-?'th relay, and ith relay-to-destination links by /o, fi, and gi, 
respectively. Suppose each link has Rayleigh fading, independent of the others. Therefore, /o, 
fi, and gi are i.i.d. complex Gaussian random variables with zero-mean and variances cTq, cr|., 
and cr^., respectively. As in [13], all links are assumed to be long-term quasi-static wherein all 
HARQ rounds of a single packet experience a single channel realization. Subsequent packets 
experience independent channel realizations. Note that such an assumption, applicable in low- 
mobility environments such as indoor wireless local area networks (WLANs), clearly reveals the 
gains due to HARQ since temporal diversity is not present. 

A. Relay Selection Strategy 

In this paper, we use selection relaying, a.k.a. opportunistic relaying [19], which selects the 
best relay among N available relays. Inspired by the distributed algorithm proposed in [19], 
which uses request-to-send (RTS) and clear-to-send (CTS) signals to select the best relay, we 
propose the following selection procedure for HARQ systems: 

• In the first step, the source node broadcast its packet toward the relays and the destination. 
Thus, relays can estimate their source-to-relay channels. 

• If the destination decodes the packet correctly, the relays would not cooperate. Otherwise, 
relays exploit the NACK signal which is transmitted by the destination to estimate their 
corresponding relay-to-destination channel. 
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• The 2th relay, i = 1, . . . , has a timer Tj which its value is proportional to the inverse of 

min{|/i|2, l^ip}. 

• The relay with maximum amount of min{|/jp, \gi\'^} has a smallest Tj. Whenever the first 
relay finished its timer, it broadcasts a flag packet toward the other relays to make them 
silent and announce her as the selected relay. 

Note that the process of selecting the best relay could be also done in a centralized manner by 
the destination. This is feasible since the destination node should be aware of both the backward 
and forward channels for coherent decoding. Thus, the same channel information could be 
exploited for the purpose of relay selection. After selecting the best relay, a feedback packet 
containing the index of the best relay should be sent from the destination toward the source and 
relay nodes. 

B. Transmission Strategy 

Let s and s denote the transmitted signals from the source and the selected relay, respectively. 
As shown in Fig. [21 during the first HARQ round, the relays and destination listen to the 
source transmit block s. At the end of the transmission, the destination transmits both the 
source and relays a one-bit ACK or NACK indicating, respectively, the success or failure of 
the transmission. The NACK/ACK is assumed to be received error-free and with negligible 
delay. Then, with the procedure given above, the best relay is selected. As long as NACK is 
received after each HARQ round and the maximum number of HARQ rounds is not reached, the 
source successively transmits subsequent HARQ blocks of the same packet. As illustrated in Fig. 
m suppose the selected relay decodes the message after HARQ round k, while the destination has 
not yet decoded the message correctly. For all HARQ rounds / > k, the source and the selected 
relay simultaneously transmit s and s, respectively. For this non-orthogonal transmission, the 
destination can be benefited from the spatial diversity using the following methods: 

1) Space-Time Code Transmission: The Alamouti code can be used to transmit the coded 
packets, hence, no interference occurs due to the simultaneous transmissions of the source 
and relay. The effective coding rate after / HARQ rounds is R/l bps/Hz, where R is the 
spectral efficiency (in bps/Hz) of the first HARQ round. Let x and x denote the Alamouti 
code transmitted signals from the source and the selected relay, respectively. The received signal 
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y at the destination can be written as follows: 

{fnX + QrX + 71, \f I > k, 
(1) 
/ox + n, if I <k. 

where the index r refers to the index of the selected relay and n is a complex white Gaussian 
noise sample with variance Nq. 

2) Beamforming Transmission: An alternative way of simultaneous transmission of the coded 
packets is beamforming. Assuming the knowledge of channel phases of /o and g^. at the source 
and the selected relay r, respectively, the transmitted message can be recovered at the destination. 
Moreover, using beamforming, we can achieve array gain of two, comparing to space-time code 
usage, in expense of higher data exchange overhead for phase estimation at the transmitters. 
However, ACK/NACK transmission from the destination can be exploited for the phase estima- 
tion of the channels. Thus, no signaling overhead is added when beamforming technique is used. 
In this case, the received signal at the destination is given by 

{l/ol-s + lfl-rls + n, if / > /c, 
(2) 
|/o|s + n, if / < /c. 

C. Average Throughput 

Two definitions of throughput are considered. A frequently used metric for throughput analysis 
is the long-term (LT) average throughput, given by [20] 

R B_ 

" Ef=o' ^out(o 

where E{/} is the average number of HARQ rounds spent transmitting an arbitrary message and 
-fout(0 denotes the probability that the packet is incorrectly decoded at the destination after I 
HARQ rounds. In the next section, we calculate closed- form solutions for the outage probability 
terms Pout(0 ^^^^ in ©. The definition in ([3]) relies on the steady-state behavior of several 
message transmissions. During this time, the probabilities Pout(0 assumed to be constant. 
This assumption is removed by considering the delay-limited (DL) throughput, which is the 
throughput of a single packet, defined by [13] 

Gdl = Y,J [^out(/ - 1) - PouM , (4) 
1=1 
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An advantage of definition dH), which does not resort to long-term behavior, is the ability to 
track slow time variations in the channels. 

In [21], the request to an automatic repeat request (ARQ) is served by the relay closest to 
the destination, among those that have decoded the message. However, distance-dependent relay 
selection does not consider the fading effect of wireless networks and leads to a maximum 
diversity of two. Therefore, in this work, the request to an ARQ is served by the relay with the 
best instantaneous channel conditions. Similar to [19], we choose the relay with the maximum 
of min{7/-,7g.}, z = 1, . . . , iV, as the best relay, where 7/. = and 7g. = \gi\^. We define 

7max = min{7/^,7gj 

= max {min {7/, , 7^ J , . . . , min {7/^ ,73^}} (5) 

where 

r = arg max {min {7^^ , 73 J } • (6) 

1=1,. ..,N 

III. Performance Analysis 

In this section, we calculate the outage probability of the HARQ relay selection system 
proposed in the previous section. Besides achieving a performance metric, outage probability 
expression is needed in both throughput definitions in ([3]) and (|4]). 

Let X denote the earliest HARQ round after which the relay stops listening to the current 

message. The outage probability for the relay channel after I HARQ rounds is given by [17] 

i-i 

Pout{l)=Y,Pout{l\l>k)Pr[x = k] 

k=l 

L 

+ $^Pout(/|/< A;)Pr[x = A;]. (7) 

k=l 

To compute Pr[x = k], the mutual information between source and relay for each HARQ round 
is given by 

Ifr = log2 (1 + ^7/.) , (8) 

where P is the average transmit power from the source and 7/^ is an exponentially distributed 
random variable with mean aj . For A; = l,...,/ — 1, x = A;if the message is successfully 
decoded by the relay at the A;th HARQ round, and we have 
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Pr[x = k]= Vr[{k - 1)1 < R, k If,. > R] 

= Fr[{k - 1)1 <R]- Fr[klf^ < R] 

= Hlfr < /Ufc-i] - Pr[7/,. < /ifc], (9) 

where 

^^. = f (2^/^= - 1) . (10) 

For k = I, . . . , L, X = k if the relay did not decode the message successfully after (/ — 1) 
HARQ rounds, and thus, we have 

Pr[x = k] = Pr[(/ - 1)/;^ <R]= Pr[7;^ < fii_^]. (11) 

From ^ and (fTTI) . Pr[x = k] can be calculated as 

p r ,1 / Pr[7/. < /ifc-i] - Pr[7/,. < /ifc], '^^ k < I, 
Pr[x = = < (12) 
[ Pr[7^, </i^_i], ifA;>/. 

A. Exact Outage Probability 

Since the index r given in ^ is dependent on channels, 7/,, and 7^,, are not independent for 
> 1. Thus, obtaining a closed-form for PDF is not straightforward. As it is seen from (fT2l) . 
for computing Pr[x = k], the CDF of random variable 7/^ is required. In the following, the CDF 
of the random variable 7/^ is derived. 

Proposition 1: Let 7/^- and 7^-, i = 1, . . . , N, be independent exponential random variables 
with means cr| and a^^, respectively. The CDF and PDF of 7/^, where r is defined as are 
given by (fT3l) and (fT4l) . respectively. 



(13) 



i=l j=l 9j ''^ i=l 



— \u p 

j=l ii i=l 



+ ^r^e \ e "'^Xlil-e )d^\ , 7 > 0. (14) 
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Proof: The proof is given in Appendix 1. ■ 
It is noteworthy that the integrals in (fT3l) and (fT?)) can be easily calculated, as all terms in the 

expansion of integrands are of the exponential form. 

Corollary 1: When cr|^ = a) and o-|^ = a| for alH = 1, 2, A^, CDF and PDF of 7/,, are 

respectively simplified as follows: 

F,j7) = (l-e''^'*")" 

-^-^e 1-e ^ ^ (15) 

f7/.(7) = -2e ^ Ml-e ^ ) 
^/ 



,e ^/S 1-e ^/ (16) 

where B{x\ a, h) = t"~^(l — t)^~^dt is the incomplete beta function [22]. 

Corollary 2: For high values of SNR, i.e. when 7 = /i^ — )■ 0, the closed form solution for 
(fT3l) can be obtained as 

j=l Ji 9i 2=1 /i ys 

From (fT3l) . Pr[x = A;] in (fT2l) can be written as 

F^, (uk~i) — F^f (ma;), if k <L 

[ F7/,.(/^«-i)> if /;;>/. 

Next, the conditional probabilities Pout(^ |/ > A;) and Pout(^ \l < k) in ^ will be calculated. 
After correct decoding of the source packet at the relay, the relay helps the source by simultaneous 
transmission according to the Alamouti code. Hence, assuming the relay transmits the same power 
P as the source, the mutual information of the effective channel is given by 

/.M = log2 (1 + ^7/0 + ^7..) . (19) 

Let Itot,k,i denote the total mutual information accumulated at the destination after / HARQ 
rounds and when x = k. For k < I, the relay listens for k HARQ rounds and transmits the 
message simultaneously with the source using the Alamouti code for the remaining (/ — A;) 
HARQ rounds. For k > I, the relay does not help the source during the / HARQ rounds. Hence, 



Submitted 11.2011 to IEEE Transactions on Vehicular Technology 



10 



'tot.fcj 



(20) 



k If^ + {l-k) Is,r,d, if = 1, 1, 
llfg, ifk = l,...,L, 

where //(, is the mutual information between the source and destination at each HARQ round 
and can be written as //^ = logg ^1 + ;^7/oj- 
Therefore, for A; > /, we have 



Pout{l \l<k) = Pr[//j„ < i?] = 1 - exp 



a 



fo 



From (|20|) . the conditional probability Pout(^ |/ > A;) can be calculated as 

Pout(/ \l>k)= Pr[/,ot,fe,i < R] 



Pr< log2 



P 



P 



P 



l-k 



< R] 



(21) 



2R/{i-k) 



Nn ~^ Nn ^/o 



k/{l-k) 



7/0 



No 
P 



7/0=0 -'79r=0 ^/u 



where /3(7/„ 



(22) 



2^/(' '-) jVo^^^^^ ~7/o ~ "p^- Due to symmetry, the PDF of random variable 7^,^, 
i.e., f^^^(7), is same as the PDF of random variable 7/^, with perhaps different mean. Thus, the 
PDF of 7g^ can be found by the derivation of Pr{7g,. < 7} in (1551) as 



N 



1 1 A' 



1 -(t^ + -^)7 



,=1 -^h 



n(i-e'^ ^ ) 



1=1 



+ 



7 _ /3 



e I e 



e ,7^0. 



(23) 



i=l 



By substituting f^^^ (7) from (dll) into dill), Pout(^ |/ > A;) is obtained. Therefore, using ([HI), (|2T1) . 
and (|22)) . the outage probability in the /th stage of HARQ process can be achieved as 

Pout(/) = 5^(F,^^(/ifc-i) - F,^^(/x,))T(/, A:) + 5^F,^^(/i,_i) ( 1 - e^fo | , (24) 
fc=i 

where T(/, A;) is defined in (l22l) . 



fc=/ 
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B. Approximate Outage Probability 

In the previous subsection, we were able to derive the outage probability in the /-th round 
of HARQ. However, a triple integral should be solved to get T(/, k), and also F^^ is in 

an integral form. To give an insight on the diversity and find ways of optimizing the system, in 
this subsection, we try to find a simpler solution for the outage probability. In the following, an 
approximation of the CDF of the random variable 7/-^ is derived. 

Proposition 2: Let 7^. and 7^-, i = 1,. . . ,N, be set of independence exponential random 
variables with mean cr|. = cr|. = af. The cumulative density function of 7^,,, where r is defined 
as can be approximated as 



F7/.(7) 



N 

n 

i=l 



1-e 



27 



(25) 



Proof: The proof is given in Appendix II. ■ 
In Fig. [3l we have compared the approximated PDF of 7/^, which is obtained by the derivation 
of CDF in (I25]) . with the simulated PDF of 7/-,,. As it can be seen from Fig. |3l for the case of a 
single-relay network {N = 1), the analytical and simulated results have the same performance. 
This is because of the fact that the independence assumption for 7/. and 7^^ becomes valid for 
= 1, and the approximation in (1251) turns into equality. For the opportunistic relaying case, 
i.e., > 1, it can be seen that the analytical curves appropriately approximate the simulation 
result. 

From (|25l) . Pr[x = k] in (fT2)) can be approximated as 



Pr[x = k] 



\ 



N 

n 

i=l 



2Mfc 



N 

n 

i=l 



12— 



for k < I, and 



Pr[x = A;] ^ 1 



N 

n 

i=l 



2M;_i 



1 - e 



n^{k), 



(26) 



(27) 



for k>l. 

Due to symmetry, the PDF of random variable 7^^, i.e., fy^^ (7), is same as the PDF of random 
variable 7^^, with perhaps different mean. Thus, the closed-form approximation for PDF of 7^^ 
can be found by the derivation of F^^^ (7) in ( |25] ) as 
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KM 



\ 



i=l 



By substituting f7g^(7) from (|i8l) into (|22)) . Pout(^ |^ > ^) can be approximated as 



Pout(/ \l>k)^ / ° ^ f,,^ (7,J ci7/o di,^ ^ r,{l, k). (29) 

Therefore, using (fT2)) . (f2T)) . and (|29|) . the outage probability in the /th stage of HARQ process 
can be achieved as 



i-i 



-1^1 



Pout(0~5Z^i(^)T2(/,A:) + 5^fi2(0 l-e^^° • (30) 

fc=l k=l \ / 

C. Upper-Bound on Outage Probability 

For calculating the minimum diversity gain of HARQ wireless relay networks when selection 
strategy in ^ is used, it is enough to derive an upper-bound on the outage probability Pout(0- 

The random variable 7/,,, which is corresponding the source-relay channel of the selected 
relay, can be bounded as 

7max < 7/. < 7max, (31) 

where 7max is given in ^ and 7^^^^ is defined as 

7max = . max {7jJ . (32) 

t=l,...,N 

The CDF of 7^^^^ can be written as 

Pr{7max < 7} = Pr{7/i < 7, 7/2 < 7, • • • , IJn < 7} 

N 



i=l 



n 1 - e • (33) 



Thus, it is easy to show that the CDF of 7/,, can be bounded as 

Pr{7max < 7} < Pr{7/. < 7} < Pr{7max < 7}- (34) 

Therefore, combining (fT2l) and (l34l) . an upper-bound on Pr[x = k] will be obtained as follows 

Pr[x = k]< Pr[7max < f^k-i] - Pr[7max < /^fe], (35) 
for A; < /. From (1531) . (1331) . and (|35]) . Pr[x = k]fox k <l can be calculated as 
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i=l \ / j=l V 



^Ki{k). (36) 



For /c > /, by combining (fT2l) . (|53]) . and (l34l) . we have 

Pr[x = /:;]< Pr[7max < /Ui-i] 



/ „ 1 , 1 



n 1-e " V-^. --^J ^A,(l). (37) 



Next, Pout(^ |/ > A;) in (l22l) can be upper-bounded as 

Pout(/ 1/ > k) 



< Pr < 7max < ^ ^ k/(l-k) '^fo 



No 



k/{l-k) 'JO p 



'} + i7/o) 

/ / P7max(7max)c?7/oC?7: 



(38) 



'7/0=0 "'7max=0 "^/o 

The PDF of random variable 7maxj i-C-^ f7max(7) be found by the derivative of Pr{7inax < 7} 
in (|53] ). Thus, we have 



X7) = E 



7.axW. Z^. . 



e 



X 



n 1-e V-/. . (39) 

Therefore, by substituting Pr[x = A;] from (l36l) and (|37]) . and Pout(^ \l < A;) and Pout(^ 1^ > k) 
from (I2TI) and (|38] ). respectively, in (|7]), an upper-bound on outage probability the /th stage of 
HARQ process, i.e., Pout(0 can be achieved. 

A tractable definition of the diversity gain is [23, Eq. (1.19)] 

^ ,. log (Pout) 

Gd = - lim - — -— , (40) 

p-*oo log [p) 

where p = Thus, in the following, we investigate the asymptotic behavior and diversity order 
of Pout(/) in m. 

From (l39l) . an upper-bound for P'y^^^i'j) can be found as 
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f,„(7)<iV7'"-'nf^ + ^), (41) 

i=l \ fi 9i / 

which is a tight bound when 7 — )■ 0. Note that in high SNR scenario, the the behavior of the 
fading distribution around zero is important (see, e.g., [24]). 

Using (|4T1) and the fact that the exponential distribution is a decreasing function of 7/„, 
-Pout(^ |/ > A;) in (|38l) can be further upper-bounded as 

Pomil \l > k) 

Combining (|7]), ((2T]) . (l36l) . (|37]) . and (l43l) . a closed-form upper-bound for the outage probability 
after / HARQ round can be obtained as 

l-\ L / -M,\ 

Pout(/) < 5^Ai(A;)vl/(Z,A;) + J]A2(/) 1 - . (43) 

k=l k=l \ J 

Furthermore, using (l26l) . and (|T7l) . another closed-form approximation for Pout(0 be 
obtained as 

l-l L / -M,\ 

nut(/) ^ 5^ ^i{k)^{l, A;) + ^ fi2(/) 1 - . (44) 

fc=i fc=i V / 

Proposition 3: Assuming a HARQ system with potential relays nodes, the relay selection 
strategy based on ^ can achieve the full diversity order of A^ + 1. 

Proof: The proof is given in Appendix III. ■ 

IV. Numerical Analysis 

In this section, the performance of the proposed relay-selection HARQ system is studied 
through numerical results. We used the equal power allocation among the source and the selected 
relay. Assume the relays and the destination have the same value of noise power, and all the links 
have unit- variance Rayleigh flat fading, i.e., cr|. = (t|. = a'j^^ = 1. It is also assumed that rate R is 
normalized to 1 . We compare the transmit SNR versus outage probability performance. The 
block fading model is used, in which channel coefficients changed randomly in time to isolate 
the benefits of spatial diversity. The simulation result is averaged over 3'000'000 transmitted 
symbols (channel realization trials). 
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Fig. m confirms that the analytical results attained in Section III for the outage probability 
have an accurate performance as the simulation results. We consider the maximum number of 
HARQ rounds to be L = 5. The outage probability at the 2nd HARQ round, i.e., Pout(^ = 2), 
is compared for two different number of relays N = 2,A. One can see the approximate outage 
probability derived in (l30l) has the similar performance as the simulated curved for all values of 
SNR. In addition, the closed- form outage probability expression in (|44l) well approximates the 
simulated results, especially in medium and high SNR conditions. Furthermore, Fig. |4] shows that 
the upper-bound expression in (l43l) is a tight upper-bound. The asymptotic outage probability 
derived in (l56l) is also depicted in Fig. |4] which confirms the full-diversity order of the proposed 
scheme. 

It is straightforward to show that the outage probability for direct transmission after / HARQ 
round is [13, Eq. (7)] 



In Fig. [5l the outage probability at the I = L = 5th HARQ round of the system with different 
number of relays are considered. After selecting the best relay, Alamouti code is employed 
in the second transmission phase. Compared to the single HARQ relaying system proposed in 
[13], the proposed HARQ opportunistic relaying system with N = 2,3,4 relays outperforms 
considerably for all SNR conditions. For example, it can be seen that in outage probability of 
10~^, the system with two relays saves around 8 dB in SNR compared to the single relay HARQ 
system. Furthermore, it can be checked that the system with relays can achieve the diversity 
order of + 1. 

A delay-limited throughput where defined in @) explicitly accounts for finite delay constraints 
and associated non-zero packet outage probabilities. It can be shown that for small outage 
probabilities, this delay-limited throughput is greater than the conventional long-term average 
throughput defined in ^ In addition to finite delay constraint, represented by the maximum 
number L of HARQ rounds, higher-layer applications usually require that Pout < Pmax, where 
Pmax is a target outage probability. The total LT and DL throughput are studied in Fig. [6] subject 
to user QoS constraints, represented by outage probability target pmax and delay constraint L. In 
Fig. [6l the total LT and DL throughputs of opportunistic relying HARQ system with = 2,4 
relays are plotted as a function of SNR and compared with the direct transmission HARQ 




(45) 
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system with L = 3, pmax = 10"^, and the following linear relay geometry: (t|^ = (t|^ = ctq = 1. 
As expected, the presence of the relays significantly increases the throughput. Furthermore, in 
agreement with [13, Eq. (5)], it can be seen that the delay-limited throughput is greater than 
the long-term average throughput. An interesting observation is that as well as the diversity 
gain achieved by the opportunistic relaying HARQ system, which is previously shown in Fig. [51 
obtaining higher throughputs are possible. This behavior underscores the importance of the 
proposed system. 

V. Conclusion 

In this paper, we proposed a throughput-efficient relay selection HARQ system over Rayleigh 
fading. The throughput-delay performance of a half-duplex multi-branch relay system with 
HARQ was analyzed. A distributed relay selection scheme was introduced for HARQ multi- 
relay networks by using ACK/NACK signals transmitted by destination. We evaluated the average 
throughput and outage error probability performance and showed that the proposed technique 
significantly reduces the multiplexing loss due to the half-duplex constraint while providing 
attractive outage error probability performance. The closed-form expressions outage probability 
were derived, defined as the probability of packet failure after L HARQ rounds, in half-duplex. 
For sufficiently high SNR, we derived a simple closed-form average outage probability expression 
for a HARQ system with multiple cooperating branches. Based on the derived upper-bound 
expressions, it was shown that the proposed scheme achieves the full spatial diversity order of 

+ 1 in a non-orthogonal relay network with parallel relays. The analysis presented here 
allows quantitative evaluation of the throughput-delay performance gain of the relay selection 
channel compared to direct transmission. The numerical results confirmed that the proposed 
schemes can bring diversity and multiplexing gains in the wireless relay networks. 

Appendix I 
Proof of Proposition 1 

First, we define the auxalary random variables rrii = min{7/-, 73.} for i = 1,2, A^. Since 
{lfi^lgi}f=i independent exponential random variables, mjS are also independent exponential 
random variables with the following CDF: 

F„^(x) = l-e , x^O (46) 



Submitted 11.2011 to IEEE Transactions on Vehicular Technology 17 

Also, using partitioning theorem, we have 

N 

Pr{7/. ^ 7} = $^Pr{(7/. ^ 7) H (r = j)}. (47) 

For j = 1,2, N, the summands of (l47l) can be obtained as follows 

Pr{(7/. ^7) n(r = j)} 

N 



Pr{(m,- ^ 7) n ( f|K < '^i)) n (7/, ^ 7)} 



i=l 



Pr{(mj ^ 7) n ( p|(mi < mj))} 



i=l 



N 



- Pr{(m,- ^ 7) n ( f|(mi < m,)) n (7/^ > 7)}. (48) 



By substituting (|48l) in (147]), we obtain 

Pr{7/. ^ 7} = 5^Pr{K- ^ 7) n ( f|K < ^i))} 
j=i i=i 

N N 

- Pr{(m,- < 7) n ( f|(mi < m,)) n (7/^ > 7)} 

j=l i=l 

= Pr{max(mi, m2, ttln) ^ 7} 

N N 

- Pr{K- ^ 7) n ( f|(mi < m,)) n (7/^ > 7)}. (49) 

7=1 i=l 

Since mjS are independent, the first term on the right side of (|49|) is given by 

N 

Pr{ max(mi, 7712, ...,mj\[) ^ 7} = Pr{[^(mj ^7)} 

AT AT 

=nprw^^7}=nF-»w- (50) 

1=1 i=l 

Also, the summand of the summation on the right side of (|49l ) is obtained as follows 
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N 



Pr{(mj- ^ 7) n ( p|(m, < nij)) n (7/^, > 7)} 



i=l 



N 



Pr{7/, > 7} Pr{(mj ^ 7) H ( p|K < ^j)) \lfj > l} 



i=l 



N 



Mifj > 7} Pr{(7<?. ^ 7) n ( f|K < 79.))} 



/■7 1 _ /5 Af 

/ _e ^Pr{n(m, </3)}rf/3 
^0 '^0, 



e ■''j 



'0 "9. i=l 



^ /•7 2_ ^ 



—e ^1 



Substituting from (gll) into dSO]) and dSJ), one can obtain the CDF in using dUl) to (|5T1) . 
Also, taking derivative of (fT3]) with respect to 7, results in the PDF of 7/^, given by (fl?)) . 

Appendix II 
Proof of Proposition 2 

For deriving the CDF of 7/^, we should first find the CDF of 7max, which can be written as 

Pr{7max < 7} = Pr{7i < 7, 72 < 7, • • • , 7Af < 7} (52) 

where 7j = min {7/^,79^} is again an exponential random variable (RV) with the parameter equal 
to the sum of parameters of exponential RV 7/-- and 7^., i.e., I/ct/. and respectively. 

Thus, assuming that all channel coefficients are independent of each others, we can rewrite 
(1521) as 

N 



i=l 




Pr{7max < 7} = n h - ^ "''^ ]■ (53) 



On the other hand, we have 

Pr{7max < 7} = 1 - Pr{min {7/,, 7g,.} > 7} 

= 1 - Pr{7/^ > 7, 73,, > 7} ^ 1 - Pr{7/^ > 7}Pr{79. > 7}, (54) 

where the last equality is an approximation as if 7/^,, and 7^^ are independent. For simplicity, we 
assume equidistance source-relay and relay-destination links, i.e., that cr^. = (t|. = af. Since we 
have assumed that 7/. and 7^. have the same statistics, using (l53l) and (l54l) . we have 
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Pr{7/. < 7} = Pr{79. < 7} ~ 1 



\ 



i=l 



n . (55) 



Appendix III 
Proof of Proposition 3 

From a Taylor series expansion, it can be shown that the first term in (l43l) is 0(l/p^^"^^). 
From (|43] ). and by representing the factor fi^ in terms of the SNR ratio p, the outage probability 
in high SNR can be written as 

Pout(/) < (56) 

where 

A(;)= 2T-1 \ 

JO i=i \ .n 9i / 

Hence, observing (l56l) . the diversity order defined in (l40l) is equal to + 1, which is the full 
spatial diversity for + 1 transmitting nodes. 
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Fig. 1. Wireless relay network consisting of a source, a destination, and A' relays. 
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Fig. 2. Example of HARQ protocol for relay selection system. The selected relay decodes message after HARQ round k. 
The source and selected relay simultaneously transmit si and s;, respectively, for all HARQ rounds I > k. In this figure, the 
destination decodes the message after HARQ round Af where k < M < L. 




Fig. 3. Comparison of the PDF of tiie received SNR at the selected relay r in a network with A*" relays. 
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Fig. 4. The outage probability Pout(0 curves of delay-limited HARQ networks employing opportunistic relaying with 2 and 4 
relays, when R = 1 bits/sec, L = 5 is the maximum number of rounds, and we consider HARQ round of I = 2. 
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Fig. 5. The outage probability performance of the proposed relay selection HARQ system versus transmit SNR in a 
with different number of relays, R — 1 bits/sec, L — I — 5 HARQ rounds, and CTj. = ct^. — aj^ = 1. 
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Fig. 6. The delay-limited (DL) and long-term (LT) throughputs of direct transmission and relay selection HARQ system 
versus transmit SNR for target outage probability pmax = 10^'*, L = 3 HARQ rounds, physical layer rate R = 1 bits/sec, and 
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